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ABSTRACT 

We analyze the spectral energy distributions (SEDs) and emission line properties of the red (J— Kg 
> 2) 2MASS AGN using principle component analysis (PCA). The sample includes 44 low redshift 
AGN with low or moderate obscuration (Nh < 10 23 cm" ) as indicated by X-rays and SED modeling. 
The obscuration of the AGN allows us to see weaker emission components (host galaxy emission, AGN 
scattered light) which are usually outshone by the AGN. The first four eigenvectors explain 70% of 
the variance in the sample. Eigenvector 1 (33% variance in the sample) correlates with the ratios 
of the intrinsic X-ray flux to the observed optical/IR fluxes and the F(2— 10keV)/F([OlII]) ratio. 
We suggest that it is primarily driven by the L/LiEdd ratio and strengthened by intrinsic absorption 
(both circumnuclear and galactic). Eigenvector 2 (18% of variance) correlates with optical/IR colors 
(B— Kg, B— R, J— Kg) and optical spectral type and depends on the contribution of the host galaxy 
relative to the observed AGN emission. Eigenvector 3 (12% of variance) correlates with reddening 
indicators obtained from the X-rays (hardness ratio, spectral index, Nh) and the narrow Ha/H/3 
emission line ratio. Their relation suggests a common absorber for the optical/X-rays lying outside 
the narrow- line region (NLR), possibly in a moderately inclined host galaxy. Eigenvector 4 (8% of 
variance) correlates with the degree of polarization and the broad Ho/H/3 ratio, indicating that dust 
which scatters the nuclear emission (continuum and the broad-line region emission) also reddens the 
broad-lincs. Our analysis shows that, although as suggested by unification schemes, the inclination 
dependent obscuration (circumnuclear and the host galaxy) is important in determining the AGN 
SEDs, the L/ LEdd ratio is the most important factor, followed by host galaxy emission. 
Subject headings: galaxies: active — quasars: general 



1. INTRODUCTION 

Principal Component Analysis (PCA) is a mathemat- 
ical tool used to reduce the dimensionality of a dataset 
without losing information. It transforms the old co- 
ordinate system where variables correlate among each 
other into a new coordinate system defined by a smaller 
number of uncorrelated variables called principal com- 
ponents. The first principal component accounts for the 
most variance in the data, and the successive principal 
components explain as much of the remaining variance 
in the data as possible (see e.g. an introduction to PCA 
by Francis k Wills 1999 or Murtagh & Heck 1987). 

Since QSOs show a plethora of emission line and 
continuum properties, PCA is ideal for studying their 
properties and understanding which are most important. 
PCA was first introduced to the AGN world by Boroson 
& Green (1992; hereafter BG92) who analyzed the opti- 
cal emission line and continuum properties of a complete, 
low redshift, Bright Quasar Survey (BQS) sample. It 
was found that the primary eigenvector (EV1) is anticor- 
related with FeIlA4570 strength (equivalent width and 
FeII/H/3 ratio), correlated with [OIII]A5007 strength 
(luminosity and peak) and H/3 FWHM, and anticorre- 
lated with the blue asymmetry of the H/3 line. Hence 
BG92 concluded that EV1 is driven by a physical pa- 
rameter that "... balances Fe II excitation against illumi- 
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nation of the narrow line region" . The BG92 EV1 was 
also found to correlate with the soft X-ray properties: 
luminosity and spectral index (Boroson & Green 1992; 
Corbin 1993; Laor et al. 1994, 1997, Wang et al. 1996, 
Brandt & Boiler 1998) and the correlations were found 
to be stronger than with the individual emission line pa- 
rameters (Brandt & Boiler 1998). This implied that EV1 
has a more fundamental physical meaning. Physical pa- 
rameters that have been suggested to drive EV1 include: 
the L/ LEdd ratio, orientation, and black hole spin. Since 
EV1 is correlated with [OIII] emission, thought to be 
an isotropic property, BG92 concluded that EV1 is not 
driven by orientation. Although the isotropy of the [O III] 
emission has since been called into question (Jackson & 
Browne 1990, Di Serego Alighieri et al. 1997, Tadhunter 
et al. 1998) the correlation of EV1 with the inclination 
independent [O II] A3727 emission, originating at larger 
distances from the nucleus than [O III] , strongly supports 
the EV1 orientation independence (Kuraszkiewicz et al. 
2000), at least for the radio-quiet subsample of the BQS. 
BG92 suggest that the L/LEdd ratio is the primary driver 
of EV1. 

BG92 EV1 was also found to correlate with the UV 
emission line properties: CIII] width, Silll] / CIII] ratio, 
CIV and NV strength, CIV line shift (Wills et al. 1999, 
Wills, Shang, & Yuan 2000, Sulentic et al. 2000), sug- 
gesting that EV1 is also linked to the physics of fueling 
and outflow and that higher accretion rates are linked to 
denser gas and possibly higher nuclear starburst activity 
(Wang et al. 2006). 

By combining observational data together with numer- 
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ical simulations, Marziani et al. (2001) and Sulentic et al. 
(2000) were able to reproduce the correlations between 
optical emission line properties found in EV1 and con- 
firm that EV1 is related to L/Lsdd, but suggest that, in 
their sample, orientation is also important. Their sam- 
ple includes both the BQS quasars from BG92 (face-on, 
low obscuration), and inclined/reddened sources of in- 
termediate type and Type 2 AGN that have FOS/HST 
spectra. 

A different way to approach the analysis of AGN prop- 
erties is spectral PCA first applied by Francis et al. 
(1992) to the analysis of more than 200 spectra from 
the Large Bright QSO Survey (LBQS; redshift range: 
0.2 < z < 3.4) or by Shang et al. (2003) to the spectra 
of PG QSOs. Their analysis is superseded by Yip et al. 
(2004) who analyze spectra of more than 16,000 QSOs 
from the Sloan Digital Sky Survey (SDSS) with a wider 
range of redshifts (0 < z < 5.4), inclinations and/or red- 
dening values. Their spectral PCA found the following 
principal components: mean spectrum, host galaxy com- 
ponent, UV-optical continuum slope, and correlations of 
Balmer emission lines. Although the results obtained 
from spectral PCA and emission line+SED PCA are dif- 
ficult to compare, it is interesting to invoke them here 
for comparison with our results. 

In this paper we analyze a sample of red (J— K5 
> 2) low-redshift 2MASS AGN with moderate absorp- 
tion Nh ~ 10 23 cm~ 2 . Their SEDs were analyzed in detail 
in Kuraszkiewicz et al. (2008; hereafter Paper I) and the 
emission line properties are presented in Kuraszkiewicz 
et al. (2008 in preparation; hereafter Paper II). In this 
paper we run PCA on the SED and optical emission line 
correlation matrix and present our interpretation of the 
eigenvectors in this moderately absorbed 2MASS AGN 
sample. It is important to state here that the results 
of the PCA analysis are dependent on the selection of 
the sample and chosen properties. We do not expect the 
same eigenvectors in the analysis of the 2MASS AGN 
sample as those found by BG92 in the unobscured, face- 
on PG QSO sample. The PG QSOs' first two eigenvec- 
tors were dominated by the parameters of the central en- 
gine (L/LEdd and accretion rate - Boroson 2002), while 
our 2MASS AGN sample, due to AGN obscuration, is 
well suited for studies of weaker components such as the 
host galaxy emission and intrinsic absorption. 

2. THE SAMPLE AND ITS PROPERTIES 

The sample includes 44 red (J-K s > 2) 2MASS AGN 
selected to have B— Kg > 4.3, and K$ < 13.8, and 
observed by Chandra (see Paper I for more details). 
The sources have K band luminosities comparable to 
those of quasars (Cutri et al. 2002) and lie at low red- 
shift (z < 0.37). The sample shows a wide range of 
observed JCs-band-to-X-ray(lkeV) slopes (1.1 < olkx ~ 2; 
Wilkes et al. 2002), and a broad range of observed po- 
larization fraction at R band (0% < P < 13%; Smith 
et al. 2002, 2003). As found in Paper I the sample's me- 
dian SED is redder than that of the blue optically /radio 
selected QSOs (Elvis et al. 1994) or the hard-X-ray se- 
lected AGN (Kuraszkiewicz et al. 2003), showing little 
or no big blue bump and implying that near-IR color se- 
lection isolates the reddest subset of AGN that can be 
classified optically. Modeling of the near-IR/optical col- 
ors and X-ray spectral analysis in Paper I implies mod- 



erate amounts of reddening (Nh < few x 10 22 cm -2 ). 
The more highly obscured AGN in this sample are most 
likely viewed at intermediate angles, where the central 
engine is obscured from our view by either circumnu- 
clear dust or dust lying further out in the host galaxy 
inclined to our line of sight. The obscuration allows us 
to see/study the contributions of weaker emission com- 
ponents such as host galaxy emission or AGN scattered 
light emission which are normally outshone by the bright 
AGN. The sample includes: 7 Type 1 AGN with optical 
spectra resembling NLS1 or BALQSOs, 11 Type 2s with 
low (for a Type 2) Nh ~ 10 22 cm -2 , and 26 intermediate 
type sources. Many of the 2MASS AGN studied here 
also show high optical polarization (P > 4%). 

3. PRINCIPLE COMPONENT ANALYSIS RESULTS 

In order to investigate the relative importance of var- 
ious parameters in determining the SED shapes and 
emission-line properties in the red 2MASS sources, we 
analyze the relations between the continuum parame- 
ters and emission-line parameters by running principal 
component analysis (PCA) on the correlation matrix. 
The first PCA run included all continuum properties 
analyzed in Paper I (Sections 4,5,6): optical and near- 
IR colors (B— R, B— K5, J— Kg), optical type, optical 
slope, a op t, from optical spectra, degree of polariza- 
tion, redshift, X-ray slope, intrinsic absorption Nh, X- 
ray hardness ratio, X-ray (IkeV) luminosity (all from 
Chandra spectral fitting), bolometric luminosity from 
SEDs, various IR/optical/X-ray flux or luminosity ra- 
tios 5 , and the emission-line parameters from Paper II: 
W\, FWHM, shift from systemic redshift of: H/3 (broad 
and narrow components), Ha (broad and narrow compo- 
nents), [OIII], [Oil], optical Fell, and the following line 
ratios: [O II] /[O III], H(3 broad /[0 III] , H(3 narrow / [O III], 

Ra narrow ^narrow ) Yla broad /H(3 broad , Fe Il/H/3 fcroa<i , 

FeII/H ( S" a "' < ™, Fell/ [OIII]. A total of 56 parameters 
were analyzed. The EW, FWHM, line shifts, most line 
ratios, and optical-to-IR and IR-to-IR luminosity ra- 
tios did not contribute to the first few eigenvectors and 
running PCA with or without them gave similar results 
(i.e. nearly identical eigenvectors). To simplify the anal- 
ysis these parameters were omitted from the final PCA 
run. The final set of parameters is shown in the first col- 
umn of Table 1. Columns 2-5 show their relation with 
the first four eigenvectors (bold font denoting parameters 
that dominate the eigenvector) and quote in percentages 
how much variance they explain. Together the first four 
eigenvectors explain 70.3% of the variance in the sample. 
Their interpretation is presented below. 

3.1. Eigenvector 1 

Most of the variance (32.9%) in the sample is due to 
the correlations between the ratios of the intrinsic 1 keV 
X-ray flux to the observed optical/IR flux: F(lkeV)/FB, 
F(lkeV)jF R , F{lkeV)/Fi, F(lkeV)/Fj, F(lkeV)/F K 
(see Fig.[T]). We choose to use the observed optical and IR 

5 Flux/luminosity ratios included: F(lkeV)/F B , F(lkeV)/F fl , 
F(lkeV)/F 7 , F(lkeV)/FK, L(0.2-0.4^m)/L(0.4-0.8Atm), 
L(0.2-0.4/xm)/L(0.8-1.6/im), L(0.4-0.8Atm)/L(0.8-l. 6^111), 

L(l-10Mm)/L(10-100Aim), L(3-60)Mm)/L(60-100Atm), L(0.8- 
1.6Mm)/L(l-100Atm), L(0.8-1.6/im)/L(60-100/W), L(0.8- 
1.6Mm)/L(3-60Atm) 
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fluxes since, in this moderately obscured/inclined sam- 
ple, these fluxes depend on intrinsic reddening and host 
galaxy emission, which are a priori unknown and hence 
difficult to subtract. The X-ray fluxes, on the other 
hand, are corrected for intrinsic absorption, which was 
estimated from the Chandra X-ray spectral fitting in the 
following way: the higher (>80) count sources were fit 
with a power-law and rest-frame absorption with both 
parameters free (hereafter referred to as the "C" fits), 
the medium (~ 30 — 80) count sources, were fit with 
a fixed power-law with an X-ray photon index T = 2 
and Nh free ("B" fits). For the lowest count sources 
(<30 counts; "A" fits), fits were made using a power- 
law with fixed T = 2 and N H = 7.6 x 10 21 cm" 2 (the 
median absorption from the C fits). The importance of 
the optical/X-ray flux ratio can be understood, for exam- 
ple, in the disk-wind model paradigm (Murray & Chiang 
f990, Proga et al. 2005, and references therein), which 
predict coupling between processes producing X-ray and 
optical/UV photons. EV f is also strongly anticorrelated 
with the F{[OIII])/F{2 - \0keV) ratio (see Fig.©. 

To understand the meaning of EVf we need to take a 
closer look at the objects with extreme values of EV1. 
Their SEDs and optical spectra are presented in Fig- 
ures [3] and [4j and their properties are summarized in 
Table 2. We find that the most negative EV1 sources 
(Fig. [3]) with lowest F(lkeV)/ 'Fb,r,i,j.k and highest 
F([OHI])/F(2-f0keV), are Type 1-1.5 with blue optical 
spectra dominated by strong Fe II and weak (relatively to 
H/3) [OIII] emission, that resemble the spectra of NLSls 
and BALQSOs. Although their low S/N X-ray spectra 
(.4 fits) do not yield an estimate for Nh, modeling of 
B— R and J— Kg colors in Paper I (Section 5) gave low 
Nh ~ 5 x 10 21 cm -2 (hence their X-ray fluxes are likely 
overestimated since Nh = 7.6 x I0 21 cm~ 2 was assumed 
for these A fits during spectral fitting) and shows that 
the SEDs at optical wavelengths are dominated by AGN 
emission (>94% of total observed flux at R band is due 
to the AGN in these sources). 

The most positive EV1 sources (Fig. |4]) with highest 
F{lkeV)/F B ,Rj,j,K and lowest F([0 III])/F(2-I0keV), 
are Type 1.9-2, with red spectra showing strong [OIII] 
and [O II] emission, and weak Fe II and H/3 emission and 
strong galactic features. X-ray spectral fitting gives Nh 
~ 10 22 cm' 2 consistent with reddening obtained from 
the modeling of B— R and J— Kg colors in Paper I (Sec- 
tion 5), that also shows high host galaxy contribution at 
optical wavelengths (>61% of the total observed flux at 
R band is due to host galaxy). 

The luminosity of the [O III] A5007 emission line, origi- 
nating from the narrow- line region (NLR) , has been sug- 
gested as an indicator of the intrinsic nuclear luminosity 
of the AGN due to the similarity of the [O III] -to- hard- 
X-ray flux ratio between Seyfert 1 and Seyfert 2 galax- 
ies (Mulchaey et al. 1994, Alonso-Herrero et al. 1997, 
Turner et al. 1997). As presented in Fig. [2] the most 
negative EV1 sources do not follow the [OIII] to hard- 
X-ray relation showing higher F([OlII])/F(2— lOkeV) ra- 
tios. This indicates that in these sources the intrinsic 
X-rays are weak. The most positive EV1 sources, on the 
other hand, have lower F([OlII])/F(2— lOkeV) indicating 
strong intrinsic X-rays. 

The differences in X-ray output relative to optical/UV 



may be explained for example by an accretion disk and 
corona model of Witt, Czerny, & Zycki (1997). In this 
model both the accretion disk and corona accrete and 
generate energy through viscosity, and division of the 
flow into optically thin (corona) and optically thick (disk) 
regions results from the cooling instability discussed by 
Krolik, McKee and Tarter (1981). The model is defined 
by 3 parameters: the mass of the central black hole, 
the ratio of the luminosity to the Eddington luminos- 
ity L/LEdd (or the accretion rate M), and the viscosity 
parameter (a V i S ) assumed to be the same in both the disk 
and the corona. The model predicts a systematic change 
in the opt/UV/X-ray SED as a function of L/LEdd- For 
high L/LEdd, the big blue bump is stronger both in lu- 
minosity and relative to the X-rays so the F(lkeV)/FB 
ratio is low (as in low EV1 sources). For low L/LEdd 
the big blue bump is weaker and the F(lkeV)/ Fb ratio 
becomes high (as in high EV1 sources; see Fig. [5j . Since 
the accretion disk and corona model predictions match 
those with EV1, they suggest that EV1 is fundamentally 
related (i.e. anti-correlated) to the L/LEdd ratio. This 
is further strengthened by the fact that the most nega- 
tive EV1 sources (with lowest F(\keV) / Fb) show optical 
spectra that characterize NLSls and BALQSOs (strong 
Fe II and weak [O III] emission) that are thought to have 
high L/Lsdd ( e -g- Pounds, Done & Osborne 1995, Boiler, 
Brandt, Fink 1996, Kuraszkiewicz et al. 2000). 

From Table 1 and Fig. Q}j,d we see that the correla- 
tion between EV1 and the F(lkeV)/Fs ratio is stronger 
than the correlation between EV1 and F(lkeV) /Fk- The 
optical/IR fluxes used to calculate these flux ratios are 
observed, not intrinsic, so they are affected by dust red- 
dening (and to a lesser extent also by host galaxy emis- 
sion) . Fluxes at optical wavelengths will be more affected 
than those at near-IR wavelengths. Since EV1 shows 
stronger correlation with the F(lkeV) / Fb ratio, which 
depends both on the L/Lsdd and reddening, and shows 
more scatter with the F{lkeV) / Fk ratio, which depends 
mostly on the L/LEdd, and less on reddening, we con- 
clude that EV1 is dominated by L/LEdd- However, dust 
obscuration also contributes as it stretches the underly- 
ing L / LEdd driven correlation, at optical wavelengths, in 
these moderately obscured 2MASS sources. If dust ob- 
scuration is inclination dependent then this result is con- 
sistent with Marziani et al. (2001) who suggest that their 
EV1 found for a heterogeneous sample of mixed type (1, 
2 and intermediate) AGN, is related to the L/LEdd ratio 
convolved with orientation. 

3.2. Eigenvector 2 

Eigenvector 2 (hereafter EV2), which explains 17.8% 
of the variance in the sample, is dominated by the ob- 
served optical/IR colors (in order of importance): B— Kg, 
B— R, J— Kg, and optical spectral type (see Table 1). In 
Paper I we found that the optical B— R color, and to 
a lesser extent the near-IR J— Kg color in the 2MASS 
AGN sample, are affected by reddening and host galaxy 
emission (and in a few, highly polarized objects, also 
by scattered AGN light). We recall the results of the 
optical/IR color modeling from Paper I in Fig. [6l The 
colors of a pure reddened AGN (median QSO SED from 
Elvis et al. 1994 reddened by Milky Way dust; long-dash- 
short-dash line) are modified by host galaxy emission 
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(Buzzoni 2005 templates), where a range of host galaxy 
strengths relative to the intrinsic, unreddened AGN at 
R band was assumed (intrinsic AGN/host galaxy ratio 
= 5 — 40). For Ay < 3 the optical/near-LR colors are dom- 
inated by the AGN. At higher Ay the B— R colors be- 
come bluer and more consistent with the host galaxy col- 
ors (represented by thick solid lines). The J— Kg colors 
are much less dependent on reddening and host galaxy, 
only visible if their contribution is strong i.e. Ay > 10 
and AGN/host galaxy ratio low (compare for example 
(Sd05;5) and (Sd05;40) curves that have a high and low 
host galaxy contribution respectively). 

In Table 2 we summarize the properties of extreme 
EV 2 objects. We find that the most positive EV2 objects 
are the bluest in B— R, implying a strong big blue bump. 
They also have the bluest B— Kg colors in sample. They 
are Type Is, with optical spectra showing strong Fe II and 
weak [OIII] emission, resembling spectra of NLS1 and 
BALQSO. Column densities derived from X-ray spectral 
fitting are low (<1.5xl0 21 cm -2 ), and the optical/IR 
colors are modeled with a pure, unreddened AGN with 
no host galaxy contribution. 

The most negative EV2 objects are red in B— Ks, B— R 
and J— Kg color. The red optical/IR colors are modeled 
with an AGN reddened by dust with Ay—Vd-W mag. 
Such reddening is sufficient to fully obscure the AGN 
at optical wavelengths, giving way to a large (98% at 
R band) contribution from the host galaxy. At near-IR 
wavelengths the reddened AGN colors still dominate in 
these sources. The red J— Kg > 3 color can only be 
produced by a reddened AGN and not the host galaxy 
which has bluer colors: J— Kg = 0.8—1 (see Fig. [6]). 

In Fig. [7] we divide the sources according to the amount 
of host galaxy contribution relative to the total emission 
at R band obtained from the optical/IR color modeling 
(numbers are from Table 7 in Paper I). We find that 
the amount of host galaxy contribution relative to the 
reddened AGN decreases with increasing EV2 values. 
This implies that EV2 depends on the amount of host 
galaxy contribution, relative to the observed/reddened 
AGN, where at one (negative) end of EV2 reside pure, 
unreddened AGN dominating the output at optical/near- 
IR wavelengths, and on the other (positive) end reside 
reddened AGN whose output at optical wavelengths is 
dominated by host galaxy emission. A stronger correla- 
tion of EV2 with the B— R color than the J— Kg color 
demonstrates that the B— R color is much more sensitive 
to host galaxy contamination than the J— Kg color, as 
shown by the color-color modeling (see Fig. [6] or detailed 
discussion in Paper I). Yip et al. (2004), who analyze 
the SDSS QSO spectra (sample with a wide range of 
redshifts, inclination angles and reddening values) using 
spectral PCA, also find that their second principle com- 
ponent is related to the host galaxy. 

3.3. Eigenvector 3 

Eigenvector 3 (hereafter EV3; 11.6% of the variance) 
is dominated (in order of importance) by: X-ray hard- 
ness ratio, X-ray spectral index Tx, column density Nh 
(measured only in the high and medium S/N Chandra 
spectra i.e. in the C and B fits) and the narrow Ha/H/3 
emission line ratio. In Fig. [8] we show two representa- 
tive correlations between EV3 and the X-ray hardness 
ratio and the narrow Ha/H/3 emission line ratio. All of 



these parameters are related to the amount of obscura- 
tion in X-rays or optical. Their relation implies a com- 
mon absorber for the optical and X-rays, lying outside 
the NLR, possibly in a host galaxy inclined to our line 
of sight (as in Polletta et al. 2008, Deo et al. 2007). This 
finding is supported by the HST/WFPC2 I-band images 
of the 2MASS AGN, which show that these sources re- 
side mostly in host galaxies with inclination angles of 
i = 50 - 75° (mean i=67°; Marble et al. 2003). We 
conclude that EV 3 is likely dominated by host galaxy 
absorption. 

The presence of redshift as a component of EV3 (Ta- 
ble 1) is secondary to the presence of Nh, for which mea- 
surements become less sensitive as redshift increases and 
the absorption is shifted out of the soft X-ray band. The 
correlation coefficient reduces significantly if Nh is cor- 
rected for redshift (assuming it is intrinsic). 

3.4. Eigenvector 4 

Eigenvector 4 (hereafter EV4; 8.0% of the variance) is 
dominated by the degree of optical polarization, broad 
Ha/H/3 emission line ratio and redshift. The relation 
between the polarization and the broad Ha/H/3 ratio 
implies that the same dust scatters the intrinsic contin- 
uum+broad lines and reddens the broad emission lines. 
Since there is no dependence on the narrow Ha/H/3 ratio 
in this eigenvector, this dust must lie outside the BLR 
but not outside the NLR, i.e., cither on the BLR/NLR 
border, or within the NLR. This is consistent with the 
conclusions reached by Smith et al. (2003) based on spec- 
tropolarimetry of the 2MASS AGN. The correlation be- 
tween EV4 and redshift is due to a correlation between 
redshift and polarization (chance probability ~ 2%; there 
is no correlation with the broad Ha/H/3 ratio - chance 
probability >60%). This may be a selection effect as 
higher redshift sources will have a higher measured de- 
gree of polarization due to the wavelength dependence 
of dust scattering and lower host galaxy dilution as one 
moves into restframe UV where the stellar contribution 
decreases. 

4. CONCLUSIONS 

PCA analysis of the SED and emission line parameters 
of a sample of 44 red J— Kg > 2, moderately obscured 
(N H < 10 23 cm" 2 ) 2MASS AGN shows that 70% of the 
variance in the sample is explained by four eigenvectors, 
each having a physical explanation: 

• Eigenvector 1 (33% of variance in sample) corre- 
lates with the intrinsic X-ray to observed optical, 
IR and [O III] flux ratios, and is strongly related to 
the L/LiEdd ratio, as shown by disk-corona models, 
and strengthed by intrinsic reddening. 

• Eigenvector 2 (18% of variance) correlates with op- 
tical/IR colors and optical type, and, as shown 
by the optical/IR color modeling, is related to 
the relative strength of the host galaxy to the ob- 
served/reddened AGN. 

• Eigenvector 3 (12% of variance) correlates with red- 
dening indicators obtained from X-rays and narrow 
emission lines, indicating a common absorber that 
affects both the X-rays and narrow emission lines 
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and probably lies in a moderately inclined host 
galaxy. 

• Eigenvector 4 (8% of variance) correlates with the 
degree of polarization and broad Ha/H/3 ratio im- 
plying the presence of dust (different from the host 
galaxy absorber in EV3) that scatters both the 
BLR and the continuum emission and reddens the 
broad lines. 

PCA analysis of the SED and/or emission line prop- 
erties of unobscured (BG92) and moderately obscured 
AGN (Marziani et al. 2001, red 2MASS sample presented 
here) points to the L/LEdd ratio as the most impor- 
tant parameter in determining the SED and emission 
line properties in AGN. Since the red 2MASS sample 
has moderate obscuration, EV1 is also slightly effected 
by reddening, as is the EV1 found by Marziani et al. 
(2001) in their heterogeneous sample of moderately ob- 
scured AGN. Successive eigenvectors pick up components 
depending on the sample chosen: unobscured AGN sam- 
ples (BG92) probe the properties of the bright central 
engine (BG92 EV2 is related to accretion rate - see Boro- 
son 2002). Moderately obscured samples, such as the red 
2MASS AGN, pick up weaker components, like the host 
galaxy emission, which become important as the redden- 
ing obscures the AGN emission. Spectral PCA of the 
moderately obscured SDSS QSOs (Yip et al. 2004) also 
finds its EV2 to be related to the host galaxy. Succes- 



sive eigenvectors in the red 2MASS sample are related 
to dust obscuration, separated into circumnuclear and 
host galaxy absorption. Our analysis is consistent with 
the orientation-dependent unification scheme for AGN, 
where a range of inclination dependent obscuration and 
host galaxy contribution (as suggested by Haas et al. 
2008) and most importantly L/LEdd ratios (as suggested 
by our PCA) can explain the properties of an AGN. 
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Fig. 1. — a, b — Correlations between the intrinsic 1 keV X-ray flux to observed B flux ratio and the intrinsic 1 keV X-ray flux to 
observed R and K flux ratios. These are examples of correlations found to compose EV1. The X-ray flux is corrected for Galactic and 
intrinsic extinction while the B,R and K fluxes are corrected only for Galactic extinction. Crosses represent the blue optically/radio selected 
quasars from Elvis et al. (1994). The 2MASS objects in all figures are delineated by their Chandra S/N classifications, with circles, squares 
and filled triangles representing sources with C (highest S/N), B (medium S/N) and A (low S/N) spectral types respectively. The dashed 
line in Fig. lb shows the B-K s =4.3 color cut used to select the red 2MASS AGN. 1258+2329 lies below this cut since SuperCOSMOS B 
magnitude (which is more consistent with the SED) is used here instead of the USNO-A2 magnitude originally used for color selection, c, 
d — Correlations between Eigenvector 1 and the ratio of intrinsic IkeV X-ray flux to the observed flux in the B and K bands. 
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Fig. 2. — Correlation between the F([0 III])/F(2— lOkeV) ratio and Eigenvector 1. X-ray fluxes are corrected for Galactic and intrinsic 
reddening. The [O III] emission— line flux is corrected for NLR reddening estimated using the narrow Hot/ Hf3 ratio. Symbols indicate X-ray 
S/N as in Fig. Q] Dashed line indicates the value of log F([0 III])/F(2 - WkeV) = -1.89±0.75 (3cr uncertainty - dotted line) found by 
Mulchaey et at (1994) for Seyfert Is and Compton-thin Seyfert 2s. 
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Fig. 3. — SED and emission line properties of AGN with the lowest EV1 values. Left column — Rcstframc far-IR to X-ray SEDs. X-ray 
fluxes are corrected for Galactic and intrinsic absorption. IR-to-UV wavelengths have not been corrected for reddening or host galaxy 
contamination. SEDs show low F(lkeV)/Fg ^ ^ ratios and blue B— R colors. Right column — Observed-frame optical spectra on Fx 
(erg s _1 cm~ 2 A -1 ) versus A (A) scale. Optical spectral type, degree of optical polarization and EV1 value have been denoted on each 
spectrum. Spectra resemble those of NLSls and BALQSOs characterized by strong Fell and weak [OIII] emission. 
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TABLE 1 

Correlation of eigenvectors with SED and emission line parameters. 



Variable 


EV1 


EV2 


EV3 


EV4 




32.9% 


17.87c 


11.6% 


o no/ 

8.07o 


Tvne 


0.0709 


0.3933 


0.0914 


—0.3100 


Polarization 


—0 0358 


—0.1634 


—0.1641 


0.4851 


B — R 


0.1543 


—0.4196 


-0.0003 


0.0992 


J — K 


—0 0563 


—0.3591 


—0.3130 


0368 


B - K 


0.0678 


-0.4521 


-0.2056 


0.0547 


N H 


-0.0057 


-0.2166 


0.3696 


-0.0887 


a ov t 


-0.2108 


0.3248 


0.0724 


0.1596 


r x 


0.0468 


0.1529 


-0.3814 


-0.2999 


z 


-0.0914 


-0.0933 


-0.3250 


0.4314 


log F lkeV /F B 


0.4035 


-0.0302 


-0.0100 


0.0291 


log F lkeV /F R 


0.3999 


0.0454 


0.0172 


0.0147 


log F lkeV /F I 


0.3894 


0.0611 


0.0415 


0.0800 


log F lkeV /Fj 


0.3861 


0.1022 


0.0230 


0.0757 


log F lkeV /F K 


0.3742 


0.1686 


0.0862 


0.0822 


Hardness ratio 


-0.0520 


-0.2803 


0.4958 


-0.0517 


Ka NLR /Uf3 NLR 


-0.1120 


0.0069 


0.3347 


0.1694 


Ua BLR /Uf3 BLR 


-0.0612 


0.0165 


0.2509 


0.5241 


[OIII]/F 2 _ wkeV 


-0.3531 


-0.0104 


-0.0227 


-0.1378 



Note. — 



Parameters that dominate an eigenvector are denoted in bold font and underlined. 
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Table 2. Objects with extreme EV1 and EV2 properties 



JN anie 


EV1 2 


lnrr TP ■ , r / ZT 1 ,, 

1Q g FlkeV 1 r B 


-Q ft 


J J£ 


Type 


p 


Nfj Notes on objects /optical spectra 


.Model 


Ay 


host 








mag. 


mag. 




% 


10 22 cm- 2 




mag. 


% 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) (9) 


(10) 


(11) 


(12) 




EV1 












Extreme EV1 objects 






















most negative EV1: 








1258+2329 


-1.112 


-5.604 


0.68 


2.043 


1 


1.0 


... (A) blue spectrum with strong Fell and no narrow lines 


pure AGN 








1516+1900 


-0.634 


-5.547 


1.30 


2.168 


1 


9.3 


... (A) blue spectrum with strong Fell and weak [OIII] emission 


(Sa05;20) 


3 


26 


0234+2438 


-0.631 


-4.676 


1.18 


2.160 


1.5 


2.6 


... (A) blue spectrum with strong Fell and weak [OIII] emission 


pure AGN 


1 



















most positive EV1: 








1021+6311 


0.706 


-2.369 


1.46 


2.585 


2 


1.0 


2.4 (C) red spectrum strong [O III],[0 II], weak H/3 and Fell 


(Sd05;20) 


22 


100 


1051+3539 


0.476 


-2.691 


1.42 


2.080 


1.9 


1.2 


0.6 (C) red spectrum strong [O III],[0 II], weak H/3 and Fell 


(Sd05;5) 


3 


61 


0157+1712 


0.446 


-3.063 


1.59 


2.706 


2 


1.4 


2.8 (C) red spectrum strong [O III],[0 II], weak H/3 and Fell 


(Sa05;20) 


7 


86 




EV2 












Extreme EV2 objects 






















most positive EV2: 








1402+2631 


0.716 


-3.396 


0.64 


2.081 


1 


0.2 


0.02 (C) blue spectrum with strong Fell and weak [OIII] emission 


pure AGN 








1714+2602 


0.657 


-3.681 


0.52 


2.210 


1 


0.9 


0.15 (C) blue spectrum with strong Fell and weak [OIII] emission 


pure AGN 






















most negative EV2: 








1040+5934 


-0.419 


-4.372 


1.50 


3.020 


1.9 


1.4 


... (A) 


(Sa05;20) 


10 


98 


1307+2338 


-0.589 


-3.857 


1.73 


3.314 


2.5° 


2.5 


... (A) reddest J — K source in sample 


(Sa05;>40) 


11 


98 



N 

% 



Note. — (1) - object name, (2) - Eigenvector 1 value, (3) - ratio of intrinsic IkeV flux to the observed flux at B band, (4), (5) - optical and near-IR colors, (6) - optical spectral 
type, (7) - degree of linear polarization at R photometric band, (8) - column density measured from Chandra X-ray spectral fitting, in brackets X-ray S/N: A - low, C - high (9) - 
notes on individual objects and spectra, (10), (11), (12) - parameters obtained from optical/IR color modeling in Paper I: host galaxy type (e.g. (Sd05;5) denotes an Sd host with 
5-Gyr stars and intrinsic AGN/host galaxy flux ratio at R band =5), AGN reddening and amount of host galaxy relative to the total observed flux at R band. 

"Spectrum used to classify the optical type of this source does not extend to H/3, however Ha emission line is narrow, hence Type 2.5. 



